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Abstract 

We have obtained the exact expression of the diffusion propagator in the time- 
dependent anharmonic potential V{x,t) = ^a{t)x'^ + bliix. The underlying 
Euclidean metric of the problem allows us to obtain analytical solutions for a 
whole family of the elastic parameter a{t), exploiting the relation between the 
path integral representation of the short time propagator and the modified 
Bessel functions. We have also analyzed the conditions for the appearance of 
a non-zero flow of particles through the infinite barrier located at the origin 
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(6<0). 



I. INTRODUCTION 

The mathematical theory of stochastic processes has proven to be not only a useful but 
also a necessary tool when studying physical, chemical and biological systems under the 
effect of fluctuations []I[ . Recent theoretical and experimental studies have shown that there 
are many situations where fluctuations plays an essential role leading to new phenomena 
induced by the presence of noise. A few examples of such situations are: some problems 
related with self-organization and dissipative structures [@J^, noise induced transitions 0, 
noise induced phase transitions 0, noise sustained patterns P, stochastic resonance in 
zero-dimensional and in spatially extended systems 0J^. 

An almost natural way to describe (Markovian) stochastic processes corresponds to an 
approach introduced by Wiener based on a sum over trajectories 0, anticipating by two 
decades Feynman's work on path integrals [|1^. This approach was later applied by Onsager 



and Machlup to some Markovian nonequilibrium processes [11|. However, even some non- 
Markov processes can be also described within this framework [|r^ . In spite of the historical 
fact, path and functional integration methods have been largely studied and applied within 
the quantum mechanical realm, while its study and application in stochastic processes is 
scarce |]13|. However, it is worth here to remark that, following earlier results |l^, path 



and functional integral methods have been widely applied in the theory of dynamic critical 
phenomena, both at equilibrium and nonequilibrium phase transitions [|l^ . 

Among others, one of the problems that has not received much attention within the path 
integral description of stochastic processes is the application of space-time transformations, 
while this kind of transformations have been largely used within the realm of quantum 
mechanics |16|. Among the few studies in this regard (see Ref. [0) a recent one refers to 



a transformation relating the diffusion propagator in a time-dependent harmonic oscillator 



with the propagator for the case of free diffusion, including a whole family of possible 



analytical solutions [|T^. Also, a formal adaptation of Duru-Kleinert-like transformations 
to the stochastic case, overcoming the main disadvantage of the direct application of such 
transformations -namely that Duru-Kleinert transformations do not link different Markov 
processes- was introduced in Ref. [|T^. However, it merits to stand out the treatment of 
transformations between "different" Wiener processes done in Friendlin's book (indicated 
with (j) in Ref. ||13|), even though they are not the exact equivalent to the Duru-Kleinert 
transformations. 

In this paper we present the exact solution for the diffusion propagator in a time- 
dependent anharmonic oscillator V{x,t) = ^a{t)x'^ + blnx. This particular choice of the 
potential can be useful to model the behavior of several physical and biological systems. 
Among them, the study of neuron models (e.g. integrate and fire models 0), stochastic 
resonance in monostable nonlinear oscillators pO| and its possible application to spatially 



extended systems 0|. Also, we can consider that the logarithmic term, whithin the path 
integral scheme (or as a Boltzmann like weight) mimics a pref actor corresponding to an 
effective energy barrier. It is clear that the possibility of having exact expressions of the 
stochastic propagator in a non-symmetrical potential can be of interest. In fact, in many 
of the above mentioned applications (specially on neuron models) the potential studied in 
this work would represent a more realistic approximation to the real behavior of the system 
under study. In other problems such us brownian motors, this asymmetry is not just an 
improvement but an unavoidable ingredient of the model. 



The approach used in this work has been inspired by the solution presented in Ref. pT 
corresponding to the exact quantum mechanical propagator in a time-dependent harmonic 
potential plus a singular perturbation. In the present case, the fact that the metric of the 
underlying space is Euclidean, allows us to obtain the exact analytical expression of the 
diffusion propagator for a whole family of functional forms of the time-dependent elastic 
parameter. 

In the next section we introduce the model we are going to study and show the proce- 



dure to be followed in order to obtain the exact form of the propagator. We also discuss 
the presence of noise induced flow of particles through the infinite barrier located at the 
origin, provided that the noise amplitude is large enough for the particles to overcome the 
deterministic drift. In section |T| we show how to obtain a family of analytical solutions. In 
the last section we make a final discussion and comment on the possible applications of the 
present results. 

II. THE DIFFUSION PROPAGATOR 

In this section we will follow, and adequately adapt, the results of the paper of Khandekar 



and Lawande |2T|. Our starting point is to consider the following Langevin equation 

x = h{x,t)+m, (1) 

where ^(t) is an additive Gaussian white noise [|I[]. That is it fulfills the conditions {^{t)) = 
and {^{t)^{t')) = 2D6(t — t'). This equation describes the overdamped motion of a particle 
in a time-dependent potential. In this work, we will consider the force term h{x,t) = 
—a{t)x — b/x which, through the relation h{x,t) = —dV/dx, corresponds to the following 
potential 

V{x,t) = -a{t)x^ + b\nx. (2) 

This potential is defined for a; > and, whenever 6 < and the elastic parameter a{t) is 
positive, it corresponds to an anharmonic monostable system composed by a time- dependent 
harmonic oscillator plus a logarithmic term which is singular at the origin. 

As it will be shown later, even in this monostable situation, the noise could be able to 
induce a flow of particles through the infinite barrier located at the origin, overcoming the 
deterministic drift whenever D > —26 holds. In fact, the meaningful condition related to 
the conservation of particles inside the system (zero flux at x = 0) is D < —2b. 

In this work we will relax the monostability condition, allowing for the time-dependent 



elastic parameter to take negative values. We will show that in this extended situation an 
asymptotic probability distribution can be reached whenever the elastic term satisfies, 



lim 

to— >-00 J((j 



/ a{s)ds = oo Vt. 

Jtn 



(3) 



We will say that in this case the potential is strongly attractive. 

The path integral representation of P{xb,tb \ Xa,ta), that is the transition probability 
associated with this Langevin equation, is given by (see for instance the book indicated as 



"i" in 



ryXi,, tb I Xfl,, taj 



V[x{t)\ exp 



L{x{t),x{t), T)dT 



(4) 



Here the stochastic Lagrangian or Onsager-Machlup |]11| functional is given, in a midpoint 
discretization, by 



.X 1 r . , / m9 1 dh(x,t) 
L{x,x,t) = —[x - hix.t)]' + 2^^- 

Replacing the actual form of h{x, t) the previous expression can be expanded to yield 

d^ 



(5) 



L = L. 



dt 



(6) 



where $ corresponds to 

$(t) 
with arbitrary to? and 
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I a{T)dT H Inx 

Jto D 



a(t)x'^ 
2D 



(7) 



2D 



x^ + (a(t)2 - a(t)) x^ + {b + D) 



x^ 



(8) 



Hence the path integral in Eq. (^) adopts the form 



(9) 



with 



J\ yXfjj t^b I -^a; ^aj 



Xj, 



P[x(t)] exp 



-2zjrr+"M^^+(''+^)^ * 



(10) 



and uj(t) = a{tY — ait). As usual, the path integral in Eq. ([T0|) is defined in a discretized 
form by 

7V-1 



K{xb, U I Xa, to) = ^irn^ An ... exp\-Y^ Sj{xj, Xj^i) I JJ dxj, 

_ JV 

with Ne = ta — tb, tj = ta + je, xq = Xa] xn = Xb] An = pvrDe:] 2 ^ and 



(11) 
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(12) 



Up to first order in e (exploiting that e <^ 1) we can use the following asymptotic form 
of the modified Bessel function 



u 



exp I j - i [oiyf - i) ^ . o,.^: 



'2txu 



ie(h) 



u 



(13) 



Using the last expression with u = Xj Xj^i / D the propagator of Eq.(|TT]) may be cast to 
the following form 



K{xb,tb I Xa,ta) = lim 



1 \^/2 



N^oo \2'nDe 
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X Y\_ exp 
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The last expression can be rewritten as 



K{xb,tb I Xa,ta) = exp 
where, 



-(3 



{4 + 



Xi 



lim 3 



N 



N-l 



Y[ e "-^ ^J le {PxjXj-i) Xjdxj, (15) 



1 



«7 = /5 I 1 + y^j- 1 ; < J < iV - 1; and /5 = — . 

Now, in order to perform the integrations of Eq. (p!5|) , we can use the equality known as 
Weber formula p2|,p3|, which is given by 
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/ e °^' lei ax) leibx) x dx = — exp 
Jo 2a 



a^ + b^ 
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and is valid for 3fJ {&) > —1, 3f?(a) > (here both conditions are fulfilled). The final result 



IS 



K{Xb,tb I Xa,ta) = y/XaXt liui a^ e^^^"'"^'"''''') J^ (^aNXaXb) 



N^oo 



(16) 



where the quantities a^, Pn and g^r are defined in appendix A. These quantities are related 
to a function Q(t) that obeys the equation (as usual, we have indicated time derivatives 
with dots) 



Q{t) - u;{t)Q{t) = 0, 

with the initial condition Qo = Q{ta) = 0. In the limit e 
appendix A) 

,. 1 Q{ta) 

hm a^ = --— - — -, 



(17) 



(A^ ^ oo), we find that (see 



(18) 



lim P7V = hm - - Q2(4) / ___ 
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lim qn 

N~*oo 



1 Q{tb 



(20) 
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To calculate the second limit, it is necessary to solve Eq. (p!?!). As it was shown in Ref. 
||18|| , the complete solution of Eq. (0) can be reduced to quadratures, with the general form 
given by 



Q{t) = hntSt) + k2St^it), 



(21) 



where 



7^,^(t) = e-/*>(^^'^ 



(22) 



-/>(^)'^^ /■* 2/:aWd, 



Jta 



dr. 



Hence, the solution fulfilling the initial condition Qo = Q{ta) = is 



Q{t)=Q{QStSt). 



(23) 



After replacing this solution into the expressions for oat, pn and qN (see appendix B), 
we finally arrive to a completely analytical expression for the transition probability 
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which can be further simplified as 
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(25) 



It is straightforward to check for some particular choices of a{t) and h that the last 
expression fulfills the corresponding Fokker-Planck equation. Albeit not so simple, we have 
also proved it for the general case. The last expression also indicates that, in order to have 
the explicit form of the propagator we only need to obtain the function Q{t) (the solution 
of Eq. (|l^) given in Eq. (pT])) for the problem under study (that is, for a given form of the 
function a{t)). We will provide a family of solutions for a rather general form of the function 
a{t) in a subsequent section. Before, we will discuss the possibility of finding a net current 
at the origin. 

A. Flow through the infinite barrier 

Let us first evaluate from Eq. (|25|) the asymptotic probability distribution, that is 



P{x,t) = lim P{x,t\ Xa,ta) 



(26) 



In the strongly attractive case (Eq. (^), it can be easily shown that St^(t) diverges 
and that TZtai^) goes to zero as t^ -^ — oo. Thus, we will make use of the expansion of the 
modified Bessel function for small argument [E2|, 



Replacing this expansion into Eq. (125|) we get 



^^W = F77I^(J) +Oiz'^'). (27) 



P{x,t) 
where g{t) is defined as 



lim [x,7^i„(t)]I^+5- 

ta^ — OO 



.-4^+h+e / ^2 



exp -TTTT^TV ' (28) 



Tie + l)i2Dgit)y+^ ' V 2D^(t), 



^(t)= hm 5,„(t)7^,„(t). (29) 



It is clear that unless the condition 

1/2 + b/D + e = (30) 

holds, the system cannot reach an asymptotic probability distribution. In fact, the term 
between the square brackets in Eq. ( pSf ) depends on the initial condition. Furthermore, it 
can be shown that the normalization of Eq. (^) gives a vanishing function of t unless the 
previous condition holds. Note that Eq. ( PPD implies, through Eq. (|T^, D < —26. This 
condition gives the maximum value of noise amplitude for the particles to be confined inside 
the interval (0, cxd]. 

This encourage us to show explicitly the existence of a noise induced probability current 
through the infinite barrier at the origin when D > —2b. Before giving a rigorous deduction, 
let us state a simple argument which provides some clue about the underlying physical 
mechanism governing this flow. As it is clear from the Langevin equation, the particle 
is subjected to both deterministic and stochastic forces. If we analyze separately both 
contributions to the particle movement near the origin, we obtain for the deterministic 



trajectory Xd{t) = \J —2ht. Comparing this result with the well known diffusive behavior, 
where the uncertainty on the particle's position grows as x^ = vDt, we reobtain the previous 
condition D > —26 for the possible appearance of noise induced leakage of particles. 

9 



The probability current at the origin J(x = 0,t | Xa,ta) can be evaluated from the 
associated Fokker-Planck equation. In the case D > —26 we obtain: 



/ b D d\ 

V X 2 ox 



x=0 \ / 



(31) 



where it can be shown that Jxa,ta{t) is a positive function of time for any given initial 
condition. Therefore, we have obtained a non-zero negative current, as previously statedQ. 

B. Asymptotic probability distribution 

In the case D < —2b there is no probability leakage. In fact, the asymptotic probability 
distribution can be obtained from Eq. (|32D and is given by: 

which can be easily shown to be normalizable. 

It is worth to study how the properties of the elastic parameter function a{t) influences 
the behavior of the function g{t), which reflects the time evolution of the width of the 
probability distribution. First, note that from the definition of g{t) given in Eq. ( PUJ ) we 
can obtain 

g{t) = -2a{t)g{t) + 1. (33) 

From this equation it can be deduced that g{t) > Vt, as must be expected for any well 
behaved probability distribution. In addition, it can be proved that in order to confine 



the particle in a small region of width Jg{t) ~ y^ an attractive force of order a{t) ~ 1/e 
is needed. On the other hand, a small attractive force of order a{t) ~ e, gives a broad 



^This somewhat counter- intuitive flux has an interesting quantum counterpart in the "fall to the 
center" effect studied by Landau p^ . 
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distribution with g(t) ~ 1/e. The hmiting cases for P{x,t) corresponding to an unbounded 
spreading (a(t) ^ ^ git) -^ oo), and to an asymptotically approach to a 6{x) distribution 
(a(t) ^ oo ^ g[t) -^ 0), can be also obtained. 

From the previous paragraph it is clear that even in the strongly attractive situation 
(see Eq. (^) the probability distribution may exhibit an unbounded spreading. In fact, we 
have already shown that even in the monostable situation a{t) > 0, but where the strength 
vanishes in time, g{t) grows indefinitely. Therefore, in order to obtain a non-divergent width 
of the probability distribution, the conditions on the attractive term have to be stronger than 
the one imposed by Eq. (|^). We may infer that the localized-probability condition should 
be related to a non- vanishing attractive strength of the time averaged potential. 

C. General localization conditions 

Let us discuss the set of conditions that ensures the asymptotic localization of the proba- 
bility distribution. From the analysis of Eq. (|33D it is clear that in order to guarantee a non 
divergent git) the elastic parameter should have the following properties. Its accumulated 
strength is positive, i. e. 

f air) dT = c>0 Vt (34) 

Jti 

where c is an arbitrary constant and ti is the nearest time which fulfills the previous equation. 
This condition is clearly fulfilled if the potential is strongly attractive. We may also infer 
that the accumulated attractive effect (Eq. (|3^)) should be non- vanishing. In other words, 
the elapsed time where the accumulated strength reaches the given constant c is bounded, 
that is 

t-ti = At < At„ Vt (35) 

where At„ = Atu(c) is the mentioned upper bound for the elapsed time. It is evident that 
condition ( PSJ ) is more restrictive than the one imposed by Eq. (^). 
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In the following section we will provide a family of examples where the probability is 
asymptotically localized. 

III. A FAMILY OF ANALYTICAL SOLUTIONS 

As already mentioned, to obtain the final expression for the diffusion propagator, we 
must first solve Eq. ( [T7| ) for a given choice of the function a(t). Because the frequency uj{t) 
depends only on the harmonic term of the potential, we can make use of any known solution 
of the simpler harmonic case. A method to generate a whole family of analytical solutions 



has been proposed in Ref. |T^ for the time-dependent harmonic oscillator. In order to reach 



such a goal the elastic parameter was written in the following form 

aW = /W + ^i|| (36) 

This allows us to find the corresponding independent solutions of Eq. ([T7|) 

_ smh{F{t)) 



,,(,) . SS^m, (38) 



where F(t) = /^*^ f(s) ds, indicating that /(t) must be an integrable function. The solution 
which satisfies the initial condition Q{ta) = reduces to 

Q(,).Q(,„)5E!i(£E=£M. (39) 



With this result, the transition probability in Eq. (pSj) adopts the analytical form 

fih) " ' DsmhiF{t,)-FiQ) 



' 1 b' 



P{xt,h I Xa,ta) = ei-^-o\^'^'^>-'''^» i^,^] Xr '^ X, 



xh{b) 



XaXb \jf{U)f{ta) 

D smh{F{U) - F{Q) 



"""^l 2Dsinh(F(t,)-F(t,)) ;• ^''"^ 

Hence, we have obtained a completely analytical expression for the propagator in Eq. 
that only depends on the choice of the elastic parameter a{t). 
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IV. FINAL REMARKS 

In this work we obtained the exact expression for the diffusion propagator in the time 
dependent anharmonic potential V{x, t) = |a(t)x^ + b\nx for a rather general choice of the 
elastic parameter. The knowledge of the exact form of the propagator can be useful to model 
different physical and biological phenomena. Particularly interesting problems, suitable to 
be studied taking advantage of this results, are realistic non-symmetric neuron membrane 
potentials ||23| and the phenomenon of stochastic resonance in a monostable zero-dimensional 
potential, in spatially extended systems, and in several neuron firing models, among others. 
A complete and recent review of these stochastic resonance topics can be found in the work 
of Gammaitoni et al 0. On the other hand, the knowledge of the exact propagator in 
the indicated time-dependent anharmonic potential can be useful as a benchmark to test 
approximate numerical or analytical procedures. Among them we can refer to some of the 
problems discussed in [^ . 



Among the several studies of stochastic resonance in monostable systems, it has been 
shown using scaling arguments and numerical experiments, that the signal to noise ratio is a 
monotonically increasing function of the noise amplitude |^. By contrast, it is quite clear 
that this increase in the response of the system cannot be unbounded. We shall present 
elsewhere our own results concerning the phenomenon of stochastic resonance in a system 
described by the potential discussed in the present paper. In this regard, the relation between 
the maximum noise amplitude and the deterministic force near the origin proved to be a 
meaningful cutoff for the increase of the response. 

It is worth to remark here that the limit 6 ^ is a (kind of) singular one. The naive point 
of view will be that, in such a limit, the form of the propagator in Eq. (^51) shall reduce to 
the one corresponding to the case of the harmonic time- dependent potential V{x) ~ |a(t)x^. 
However, this limit corresponds to a harmonic time-dependent potential for x > with an 
absorbing boundary condition at x = 0. Then, it should be possible to re-obtain the limit 
6 ^ of the diffusion propagator found in this paper (Po(a^65^b I Xa.ta)) from the one 
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obtained in Ref. |jT8[| for the harmonic time-dependent case {Ph{xb,tb \ Xa,ta)) simply as 
Po{xb,tb I Xa,ta) = Ph{xb,tb \ Xa,ta) - Ph{-Xb,tb \ Xa,ta). It Can be casily proved that this 
is indeed the case. 
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APPENDIX: A 

The quantities ajy, p^ and qjy are defined according to 

«Af = P 11 7^; 7i = "1, 7j = "i - 



/ii27/ - - '^ ^ 47,-1 



P7V= - 77 + !L XT' (In = -7; + 



2 ,^147/ '" 2 47^_i' 

A = /?; /5, =/?n^> (Al) 

fc=i 27fc 
In order to determine the limiting (when A^ -^ 00) values of a^, p^, and q^, it is usefull to 
define the following auxiliary quantities 

A, = |7. ; A,= n^, (A2) 

with aj as defined after Eq. (|T3|) and 7, adopting the form 

that allows us to obtain the following equation for A, 
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A, ^2(1 + '^]- ' (A4) 



2 / Aj_i 
If we now define that A,- = -^^, the last equation can be rewritten as 

Qj+i - 2Qj + Qj^i = Ujs'^Qj, (A5) 

which, in the hniit N -^ 00 (and e ^ 0), becomes Eq. (|T7|), with the initial condition 
Qo = Qita) = 0, that follows from Eq. (^). 

Finally, we can express the coefficients a^, p^, and q^ as functions of the new variables 

/PA /3^1 pQl ~ Qo n fQl—Qo\ 1 /.„N 

aN = PAn-1 = P-T^ = P^-^ = Pe -— (A6) 

pL ^^ Ql \ (A7) 

2 y j=i Qj+iQj j 

Pe (Qn - Qn-i\ / .qs 



APPENDIX: B 



The replacement of the general solution for Q{t) indicated in Eq. (^31) into Eqs. 



r^,|^), leads us to obtain the limiting values of a^, Pn and qn- For p^ we find 



1 -1/1 1 1 , /T. N 

lim vm = lim — — — I r? r? . (Bl) 

Making a Taylor expansion up to second order in e of the last denominator we can 
calculate the limit in Eq. ( pl|) yielding 

hm PN = ^i ^-r—T- + «(4) ■ (B2) 

The expressions for a^ and q^, in terms of the explicit form for Q{t) results in 

-1 ( e/'>«* 

1 e/*>«'^^ 
lim Oat = F? • (B4) 

'J In. 
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